Meta-material photonic funnels for sub-diffraction light compression and propagation 
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We present waveguides with photonic crystal cores, supporting energy propagation in subwave- 
length regions with a mode structure similar to that in telecom fibers. We design meta-materials for 
near-, mid-, and far-IR frequencies, and demonstrate efficient energy transfer to and from regions 
smaller than 1/25-th of the wavelength. Both positive- and negative-refractive index light transmis- 
sions are shown. Our approach, although demonstrated here in circular waveguides for some specific 
frequencies, is easily scalable from optical to IR to THz frequency ranges, and can be realized in 
a variety of waveguide geometries. Our design may be used for ultra high-density energy focusing, 
nm-resolution sensing, near-field microscopy, and high-speed all-optical computing. 



While light emission by atoms, molecules, quantum 
wells, quantum dots and other quantum objects occurs 
from regions smaller than 10~ 8 to- wide, light propaga- 
tion takes place on much larger, wavelength scales, rang- 



ing from 10 



for GHz radar radiation to 10 



for 



optical light. Such a huge scale difference, effectively in- 
troduces a fundamental barrier between one's ability to 
generate light and to subsequently guide this light. Ineffi- 
cient cross-scale coupling limits the resolution and sensi- 
tivity of near-field microscopes0|3], prevents fabrication 
of ultra-compact all-optical processing circuits^, 0], in- 
tegrated optoelectronic devices 0, @ and other photonic 
systems. Here we present a class of compact waveguides 
capable of effective radiation transfer to, from, and in 
subwavelength regions. Our approach, although demon- 
strated here in circular waveguides for some specific fre- 
quencies, is easily scalable from optical to IR to THz fre- 
quency ranges, and can be realized in a variety of waveg- 
uide geometries. 

A typical waveguide structure comprises at least two 
different regions: the transparent core that supports en- 
ergy propagation, and the cladding that is used to pre- 
vent guided light from escaping the core. Although the 
design of both core and cladding may vary among differ- 
ent waveguides H, 13, E3, EH O, EE Q E3 , the wave 
propagation in the confined spaces has some universal, 
design- independent properties. Specifically, the electro- 
magnetic radiation in any waveguide forms a series of 
system-specific waves: waveguide modes. The wavevec- 
tor of each mode in the direction of mode propagation k z 
is related to the fre que ncy through the following disper- 
sion relation |l3l Il4 . : 



k = 



(1) 



where e and v are mode-specific propagation constants. 
In waveguides with non-magnetic isotropic homogeneous 
core e is the dielectric permittivity of the core, and the 

2 2 2 

parameter v = 1 — relates the frequency u), the 

speed of light in the vacuum c, the mode confinement 
radius R, and a generally non-integer mode number m 
(l m l ^ !)■ The phase velocity of the mode is given by 
the effective index of refraction n = ±v/ii7. 
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FIG. 1: (a) Schematic geometry of a waveguide with PC 
core, (b) Unfolded dispersion diagram of the fundamental 
mode in three different photonic crystals (normalized to crys- 
tal period A). Black line: TM mode in a waveguide with 
homogeneous core (ei = £2; ai = 02); note the frequency 
cut-off at \k z A\ <C 1 corresponding to diffraction limit; this 
behavior is similar to that of TE modes in PC structures 
described in the text. Blue line: TM mode with positive re- 
fractive index in PC with ei = 12; €2 = —3.8; di = 0,-2 = A/2. 
Red line: TM mode with negative refractive index in PC 
with ei = 12; £2 = — 73; ai = a2 = A/2. Our design re- 
lies on |fc 2 A| -C 1 region (highlighted area and the inset in 
(b)), which is weakly affected by crystal disorder or other im- 
perfections. The low-frequency regime (dotted region in the 
inset) is not realizable in nature since e(cu — > 0) has to be 
positive^], (c) Confined mode propagation in PC waveg- 
uides. Absolute values of real (solid lines, upper graph) and 
imaginary (dashed lines, lower graph) parts of the wavevector 
are shown for TM01 mode in the Si-core (black), and Ag-Si 
(red), Si-SiC (blue), and InGaAs-AlInAs (green) PC struc- 
tures (see text). The wavevector and radius are normalized to 
PC period and free-space wavelength respectively. Note that 
the negative-refraction modes propagate only in sufficiently 
thin waveguides Jl6jl . 



Since the mode propagation is possible only when the 
effective refractive index is real, the product of tv is 
required to be positive for such a propagation to take 
place. For the isotropic homogeneous system, the con- 
dition n 2 > is equivalent to e > 0, v > 0;?i > 0. 
Thus, there exists a minimal critical radius of a waveg- 
uide supporting at least one confined propagating mode, 
Rq ~ 7rc/(w v / e) = A/(2y / e). The systems with radius 
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R < Rq, formally described by e > 0, ^ < 0, reflect al- 
most all incoming free-space-like radiation[l|,|23|], and are 
not suitable for energy transfer into subwavelength areas. 

The properties of the waveguide modes can be con- 
trolled by cither changing the mode structure (modify- 
ing the parameter m), or by varying the waveguide core 
material (modifying the parameter e). Since there ex- 
ists only a limited control over the dielectric permittivity 
of natural materials, the majority of modern subwave- 
length waveguides implement the former technique, and 
make use of the special type of modes at the metal- 
dielectric interfaces, known as surface waves, which for- 
mally correspond to m? < 0. Although these modes 
may in principle be used to transfer radiation to nano- 
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their spatial struc- 
ture is fundamentally different from that in free-space 
waves and telecom fiber modes. This structural differ- 
ence requires non-trivial coupling mechanisms to convert 
the radiation from free-space to surface waves, typically 
associated with substantial coupling losses [15[ Il8| . 

Here we present an alternative approach to compress 
and propagate the radiation below the free-space diffrac- 
tion limit. Instead of changing the structure of the 
modes, we propose to change the waveguide itself. We 
use a periodic array of thin dielectric (e > 0) and "metal- 
lic" (e < 0) layers, widely known as a ID photonic 
crystal (PC) medium, as a meta-material waveguide core 
(Fig.la). 

In the case of PC layers perpendicular to the direction 
of mode propagation considered here, all modes of the 
system can be separated into two fundamentally differ- 
ent groups. The modes in the first group, known as TE 
waves, have their electric (E) vector parallel to the lay- 
ers, the modes in the second group (TM waves) have their 
magnetic (H) vector parallel to the layers. Similar to the 
case of a homogeneous waveguide described above, the 
frequency and wavevector of a wave in a PC-core fiber 
can be related through the dispersion relation: 



..TM 
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cos[fc z (ai + 02)] = cos(fciai) C0s(fc2O2) — 
7sin(fciai) sin(fc2a2) 



(2) 



where a\ and 02 are thicknesses of the layers in the pho- 
tonic crystal, t\ > 0, and 62 < are their permittivities, 
k\ = e^/c 2 - n 2 m 2 /R 2 , k 2 = 6 2 uj 2 /c 2 - n 2 m 2 /R 2 , 
and the parameter 7 is equal to jte — ^(f 1 + f 2 ) 
and lTM = + Sff) for TE and TM modes 

respectively[20f. The properties of several modes in typ- 
ical PC systems are illustrated in Fig. 1. 

In the case when the period of the system is much 
smaller than the wavelength and the waveguide radius 
(|fciai| <C 1, \k 2 a 2 \ < 1 and \k z (ai + a 2 )\ < 1), Eq. © 
becomes identical to Eq. Q with polarization-specific 
propagation parameters e and v, given by 



eie 2 (ai + a 2 ) R 2 u> 2 
1 n 2 m 2 c 2 



(3) 
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In a way, the PC core plays the role of a homogeneous but 
anisotropic uniaxial meta-material with its optical axis 
parallel to the direction of mode propagation[2lJ. The 
existence of propagating modes in these systems can be 
once again related to the effective index of refraction n. 
The propagation of TE modes is completely analogous 
to the wave propagation in isotropic systems described 
earlier. In contrast to this behavior, the PC structure 
can support TM waves in two fundamentally different 
regimes. The first regime, described by eJ,Q > 0; > 
corresponds to positive index of refraction, while the 
second one, e^>M < 0;vj>n < describes a negative re- 
fraction case[lfl|23,|2a,|23,|2g, unachievable in conven- 
tional fiber- and plasmonic- waveg 

uides[lElElElEl 
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Both n > and n < structures may support wave 
propagation in highly-confined areas (Fig. 1). Indeed, the 
effective refractive index of a propagating TM mode in 
substantially thin (i? <C A) strongly anisotropic system 
is inversely proportional to a mode confinement scale R. 
The decrease in R is accompanied by a decrease of "inter- 
nal" wavelength A/|n|, virtually eliminating the diffrac- 
tion limit in proposed structures. 

Furthermore, the PC waveguides with different refrac- 
tive indices can be combined together, opening the door 
for the effective phase manipulation of light propagating 
in highly-confined areas. The possibility of such a ver- 
satile light management on nanoscale is one of the main 
points of this work. 

The self- adjustment of PC waveguide modes to the 
waveguide size, accompanied by compatibility between 
the mode structure in PC waveguides, telecom fibers, 
and free-space makes the PC systems ideal candidates 
for effective energy transfer between macroscopic wave- 
propagation systems and nano-scale objects. In these 
coupling structures, called photonic funnels, the size of 
the PC waveguides gradually varies along the direction of 
mode propagation, squeezing the light into nm-sized ar- 
eas much like a conventional funnel squeezes water into 
a thin bottleneck, which is another main point of this 
work. 

The efficiency of energy compression in photonic fun- 
nels can be related to adiabaticity parameter 5 = 

d ^ 1 ^ z ^ > [H Ej| , that defines the reflection by the fun- 
nel structure, and absorption in the system. Increase of 
the funnel's length typically reduces reflection loss, but 
increases absorption. 

We now illustrate the light propagation in PC struc- 
tures described here. Although our design does not 
impose any limitations on the waveguide geometry or 
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FIG. 2: TM01 mode propagation (A = llpim) from cylindri- 
cal Si waveguide (z < 0) to a conical one with circular core; 
(a) the wave in a system with Si core. Note the reflection 
from the point where radius reaches cut-off value _Ro — 1.2/im 
(dashed and dotted lines). Only 10 -10 of energy is transmit- 
ted from R — 2.3/im to R — 0.35/im ~ A/31. This behavior is 
similar to that in tips of near- field microscopesQ,[2. (b) field 
concentration in Si-SiC PC funnel described in the text: 13% 
of energy is transmitted to 7? = 0.35/im, 16% is reflected back 
to Si (z < 0) waveguide. Note the dependence of the internal 
wavelength on the radius. 



waveguide boundary material, here we use conical waveg- 
uides with perfectly conducting metallic walls (see |29j) 
and circular cores, having adiabaticity parameters 6 ~ 
0.1 . . . 0.3, and defer the optimization of a photonic fun- 
nel geometry to future work. To demonstrate scalability 
of our design we calculate light propagation through the 
PC formed by (i) 100-nm-thick layers of SiC and Si with 
operating wavelength A = ll//m, (ii) 15-nm-thick layers 
of Ag and Si with operating wavelength A = 1.2/im, and 
(iii) 75 rom-thick InGaAs layers doped with electrons to 
10 19 cm -3 , and 150 nm AlInAs barriers with operating 
wavelength A = 20/im. (Fabrication of these structures is 
accessible through standard MOCVD, e-beam writing, or 
MBE techniques; see Refs. [l(j, [ljj, |3fJ, |3l| and references 
therein) . 

To compute the light propagation in conical struc- 
tures, we represent each structure by an array of cylin- 
drical segments (in a way, this approach also accounts 
for the effects related to finite roughness of waveguide 
walls, unavoidable in any experiment). The typical ra- 
dius step in our calculations is 10 _3 A. We then represent 
the field in each cylindrical segment as a series of modes 
of a circular waveguide. We use Eas. l2l3fl to calculate 
a mode propagation in each segment. In these calcu- 
lations we use experimental values of permittivities for 
Ag, Si, SiC, and AlGaAs[32 and use Drude approach to 
describe InGaAsll5j. Finally, we use the boundary con- 




FIG. 3: Negative refractive index systems, (a-b) homoge- 
neous Si waveguide (z < 0) coupled to Ag-Si PC core cone 
as described in the text; A = 1.2/im, (a) TM01 mode trans- 
fer from Si to a PC structure: 24% of energy is transmitted 
from R = 135nm to R = 40nm ~ A/26, 21% is reflected 
back to Si waveguide, (b) TM11 mode transfer from PC sys- 
tem (z > 0) to Si waveguide: 11% of energy is transferred 
from R = 75nm ~ A/16 to R = 230nm, 13% is reflected 
back, (c-d) TM01 energy transfer from AlInAs waveguide 
(z < 0) to passive (c) and active (d) InGaAs- AlInAs PC-core; 
A = 20/um. The passive structure transmits 6% of radiation 
from R — 2.4/xm to R = 0.45/im ~ A/44. The material gain 
in active system (eAiinAs — 10 — i) compensates for losses 
in the funnel structure and yields energy in R — 0.45/im to 
be 112% of the incident one. Similar to Fig. 2, the internal 
wavelength is proportional to R 



ditions to relate the modes in one segment to the modes 
in the neighboring segment, solving the problem of wave 
propagation through a photonic funnel. 

In Fig. 2 we demonstrate the perspectives of photonic 
funnels by comparing the energy propagation through 
mid-IR PC waveguide with positive refraction described 
above to the propagation through the Si-core structure 
with identical geometry j^. As expected, despite almost 
adiabatic radius compression, the energy in Si-core sys- 
tem reflects from the point corresponding to the cut- 
off radius of about 1.2/im. In contrast to this behav- 
ior, PC system effectively compresses energy, and the 
propagation in the structure with radius as small as 
0.35/im ~ A/30 is clearly seen. This PC provides a so- 
lution to the fundamental problem of coupling to the 
subwavelength domain, and allows transferring 13% of 
energy, which is 10 9 times better than its Si counterpart. 

The effective energy transfer across multiple scales 
in "negative-refraction" near- and far-IR PC systems is 
shown in Fig. 3. Our calculations suggest that Ag-Si sys- 
tem may be used to transfer more than 20% of energy to 
near-field zone. We expect that this number can be fur- 
ther optimized by reducing the wave reflection (currently 
21%). 
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The performance of the PC-based waveguides is lim- 
ited by the PC microstructure, and by material absorp- 
tion. The former introduces implicit inhomogeneity scale 
(PC period), where the "effective medium" approxima- 
tion [Eq. ijjjjl] breaks down. The spatial dispersion, asso- 
ciated with field inhomogeneities on such a scale, leads 
to the mode cut-off and prohibits the mode propagation 
when the radius of a waveguide becomes smaller than 
PC period. The appearance of such a cut-off is shown in 
Fig. lc. 

Material losses, on the other hand, lead to energy at- 
tenuation and limit the length of passive photonic funnels 
to ~ 10A which is acceptable for the majority of applica- 
tions of these systems: near-field tips, ultra-compact de- 
tectors, wires in all-optical circuits, etc. This limitation 
is not applicable to waveguides with active cores. Indeed, 
material absorption can be substantially reduced, elimi- 
nated, or even reversed by implementing a gain medium 
into e > regions of PC|3J|. We illustrate this approach 
in Fig. 3d by introducing moderate gain into AlInAs part 
of the far-IR structure, which can be realized via quan- 
tum cascade technology [30L l35| . 

Finally, we note that operating frequency of the pho- 
tonic funnels described here can be changed from optical, 
to near-IR, to far-IR, to THz domain by varying the PC 
composition and periodicity. The PC-based waveguides 
may be used in ultra-compact all-optical and electro- 
optical devices, near-field microscopy, and other appli- 
cations requiring effective sub-diffraction and cross-scale 
energy trans fer, as well as in a variety of nonlinear optical 
applications 36] in positive- and negative-index materials 
since the energy compression and corresponding enhance- 
ment of local field will result in the strong enhancement 
of nonlinear field moments. 

The authors would like to thank E. Mishchcnko for 
fruitful discussions 
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